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Abstract. We construct the solution of the Riemann problem for the shallow 
water equations with discontinuous topography. The system under consider- 
ation is non-strictly hyperbolic and does not admit a fully conservative form, 
and we establish the existence of two-parameter wave sets, rather than wave 
curves. The selection of admissible waves is particularly challenging. Our 
construction is fully explicit, and leads to formulas that can be implemented 
numerically for the approximation of the general initial-value problem. 



1. Introduction 

1.1. Shallow water equations. We consider the one-dimensional shallow water 
equations 

dth + d^{hu) = 0, 

dt{hu) + d^{h{u^ + 5^)) = -ghd^a, (1.1) 
dta = 0, 

where h denotes the height of the water from the bottom to the surface, u the 
velocity of the fluid, g the gravity constant, and a the height of the river bottom 
from a given level. Following LeFloch [TH] we supplement the first two balance 
laws for the fluid with the equation dta = corresponding to a flxed geometry. 
Adding the equation dta = allows us to view the shallow water equations (the 



first two equations in (1.1)), which form a strictly hyperbolic system of balance 
laws in nonconservative form, as a non-strictly hyperbolic system of balance laws 
with a linearly degenerate characteristic field. 

We are mainly interested in the case that a is piecewise constant 

/X ( aL, X < 0, 
a(x) = < „ 
^ ' \ aR, x>0, 

where , a^ are two distinct constants. The Riemann problem associated with 



( 1 . 1 1 is the initial- value problem corresponding to the initial conditions of 



(n,u,a)(x,0) ^ < ;, \ „ (1.2) 

^ ' ' ' ^ \ {hR,UR,aR), x>0. ^ ^ 

Since a is discontinuous, the system ( |1.1[ ) cannot be written in a fully conservative 
form, and the standard notion of weak solutions for hyperbolic systems of conser- 
vation laws does not apply. However, the equations still make sense within the 
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framework introduced in Dal Maso, LeFloch, and Murat (For a recent review 

see [iniisn]) 

1.2. DLM generalized Rankine-Hugoniot relations. Consider an elementary 



discontinuity propagating with the speed A and satisfying the equations (1.1). Ob- 
serve that the Rankine-Hugoniot relation associated with the third equation in (1.1) 
simply reads 

-A[a] = 0, (1.3) 

where [a] a_|_ — a_ denotes the jump of the bottom level function a, and a± 
denotes its left- and right-hand traces. Then, we have the following possibilities: 

(i) either the component a remains constant across the propagating disconti- 
nuity, 

(ii) or a changes its levels across the discontinuity and the discontinuity is 
stationary, i.e., the speed A vanishes. 

This observation motivates us to define the admissible elementary waves of the 



system (1.1 1. First of all, assume that the bottom level a remains constant across 
a discontinuity; then, a should be constant in a neighborhood of the discontinuity. 
Eliminating a from ( |1.1[ ), we obtain the following system of two conservation laws 

dth + d^{hu) = 0, 

d,{hu) + d,{h{u^+g^))^0, ^^'^^ 

Thus, the left- and right-hand states are related by the Rankine Hugoniot relations 
corresponding to (1.4) 

- X[h] + [hu] = 0, 

-\[hu] + [h{u'+g^)]^0, ^^'^^ 

where [h] := h+ — h+, etc. 

Second, suppose that the component a is discontinuous so that the speed van- 
ishes. Then, the solution is independent of the time variable, and it is natural to 
search for a solution obtained as the limit of a sequence of time- independent smooth 



solutions of (1.1 1. (See below 



Suppose that (x,i) i— > {h,u,a) is a smooth solution of (1.1). Then, the system 



( 1 . 1 1 can be written in the following form, as a system of conservation laws for the 



(now conservative) variables {h,u,a): 

dth + d,{hu) = o, 



dtu + d,{'^+g{h + a})^0, (1.6) 



dta = 0. 

Hence, time-independent solutions of ( |1.1[ ) satisfy 

(hu)' = 0, 
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where the dash denotes the differentiation with respect to x. Trajectories initiating 
from a given state (/iq, wq, ao) ^re given by 

hu = /iqUqi 

1*2 ul (1-8) 

— + g{h + a)= -j-+g{ho + ao). 

It follows from ( |1.8| that the trajectories of ( |1.7| can be expressed in the form 
u — u{h), a — a{h). Now, letting h —t h± and setting u± — u{h±), a± — a{h±), we 
see that the states {h±,u±,a±) satisfy the Rankine-Hugoniot relations associated 



with (1.6 1, but with zero shock speed: 

[hu] = 0, 



r- 

^ 2 



(1.9) 

+ g{h + a)] = 0, 



The above discussion leads us to define the elementary waves of interest, as 
follows. 



Definition 1.1. The admissible waves for the system (1.1) are the following ones: 



(a) the rarefaction waves, which are smooth solutions of (1.1) with constant 



component a depending only on the self- similarity variable x/t; 



(b) the shock waves which satisfy (1.5) and Lax shock inequalities and have 
constant component a; 



(c) and the stationary waves which have zero speed and satisfy (1.9) 



As will be checked later, the system ( 1.1 1 is not strictly hyperbolic, as was already 



observed in the previous work Recall that therein we studied the Riemann 

problem in a nozzle with variable cross-section and constructed all of the Riemann 
solutions. The present model is analogous, and our main purpose in the present 
paper is to demonstrate that the technique in [21] extends to the shallow water 
model and to construct the solution of the Riemann problem. The lack of strict 
hyperbolicity and the nonconservative form of the equation make the problem par- 
ticularly challenging. Some aspects of this problem are also covered by Alcrudo and 
Benkhaldoun |T] . For works on various related models including scalar conservation 
laws we refer to[3a[I3[Il[TSl[ni|ll|Hl[g. 

1.3. Results and perspectives. As we will show, waves in the same characteristic 
field may be repeated in a single Riemann solution. This happens when waves cross 
the boundary of the strictly hyperbolic regions and the order of characteristic speeds 
changes. We will also show below that the Riemann problem may not always have 
a solution. The Riemann problem may admit exactly one, or two, or up to three 
distinct solutions for different ranges of left-hand and right-hand states. Thus, 
uniqueness does not hold for the Riemann problem, as was already observed for the 
nozzle flow system. 

Each possible construction leads to a solution that depends continuously on the 
left-hand and right-hand states. This is a direct consequence of the smoothness 
of the elementary wave curves; by the implicit function theorem, the intermediate 
waves depend continuously on their left- or right-hand states as well as on the 
Riemann data. These results agree with [2T] which covered fluids in a nozzle with 
variable cross section. 
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In the present model, the curve of stationary wave is strictly convex. To find 
stationary waves, one needs to determine the roots of a nonlinear equation (see 
the function (p in (3.1)) which is convex and, therefore, can be easily computed 
numerically. The Riemann solver derived in the present paper should be useful 
in combination with numerical methods for shallow water systems developed in 
[21 m [m [TOl [TH [6] for which we refer to the lecture notes by Bouchut f5]. 



2. Background 

2.1. Shallow water equations as a non-strictly hyperbolic system. We now 

discuss the system (1.1 1 in the nonconservative variables U = {h,u,a). From (1.6) 
if follows that, for smooth solutions, (1.1 1 is equivalent to 

dth + udxh + hdxU = 0, 

dtu + gdxh + udxU + gdxa = 0, (2.1) 
dta = 0, 

which can be written in the nonconservative form 

dtU + A{U)dxU = 0, (2.2) 
where the Jacobian matrix A(U) is given by 

/u u 0^ 

A{U) - .9 u g 

\0 Oy 

The eigenvalues of A are 

\^{U):^u-^<X2{U):^u+^, \^{U):=Q, (2.3) 
and corresponding eigenvectors can be chosen as 

n{U) := {h,-^,Q)\ r^iU) := {h,^,0)\ 
rsiU) {gh, -gu, - ghf . 
We see that the first and the third characteristic fields may coincide: 

{\,{U)MU))^{UU)MU)) 

on a hypersurface in the variables {h, u, a), which can be identified as 

C+:^{{h,u,a)\ u^^}. (2.5) 
Similarly, the second and the third characteristic fields may coincide: 
(A2(f/),r2(C/)) = (A3([/),r3(t/)) 
on a hypersurface in the variables (/i, u, a), which can be identified as 

C- := {{h,u,a)\ u^-^/^}. (2.6) 
The third eigenvalue {X^^r^) is linearly degenerate, and we have 

~VAi([/) • ri(C/) = VA2(C/) • r2{U) - ^v^ff^ ^ 0, h>Q. 

Note also that the first and the second characteristic fields (Ai,ri), (A2,r2) are 
genuinely nonlinear in the open half-space {(/i,m, a)| h > 0}. 
It is convenient to set 

C = C+\JC^ = {{h,u,a)\ u^-gh^Q}, 



(2.4) 
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Figure 1. Projection of strictly hyperbolic regions in the (/i, u)-plane 



which is the hypersurface on which the system fails to be strictly hyperbolic. 
In conclusion we have established (cf. Figure [T]): 

Lemma 2.1. On the hypersurface C_|_ in the variables (h, u, a) the first and the third 
characteristic speeds coincide and, on the hypersurface C_, the second and the third 
characteristic 



's coincide. Hence, the system (1.1) is non-strictly hyperbolic 



The hypersurface C divides the phase domain into three disjoint regions, denoted 
by Ai,A2 and A3, in which the system is strictly hyperbolic. More precisely, we 
define 

Ai := {{h,u,a) e R+ x R x R+\ X2{U) > Xi{U) > XsiU)}, 

A2 := {{h,u,a) eR+xRxR+\ X^iU) > XaiU) > Xi{U)}, 

A+ := {{h,u,a) e A2\ u > 0}, (2.7) 

A2 := {{h, It, a) G A2I u < 0}, 

A3 := {{h,u,a) eR+xRxR+\ X^iU) > X2{U) > X^{U)}. 

The strict hyperbolicity domain is not connected, which makes the Riemann prob- 
lem delicate to solve. 

2.2. Wave curves. We begin by investigating some properties of the curves of 
admissible waves. 

First, consider shock curves from a given left-hand state C/o — {ho,uo,ao) con- 
sisting of all right-hand states U — {h, u, a) that can be connected to Uo by a 



shock wave. Thus, it follows from (1.6 1 that U and Uo are related by the Rankine- 
Hugoniot relations 



X[h] + [hu] = 0, 



(2.8) 
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where [h] ~ h ~ ho, etc, and A — X{Uo, U) is the shock speed. 

Fix the state Uq- A straightforward calculation from the Rankine-Hugoniot 
relations ( |2.8| shows that the restriction to the {h, u) plane of the Hugoniot set 
consists of two curves given by 



—o±^lih-ho)^{l + ^). (2.9) 



Moreover, along these two curves it holds 



\j h ho 

±. /-^ as ft, ho- 
V n-o 

Since the ith-Hugoniot curve is tangent to ri{Uo) at Uo, we conclude that the first 
Hugoniot curve associated with the first characteristic field is 

HiiUo): u:^ui{h,Uo)^uo^ ^{h-ho)^lQ^+J^, h>0, (2.10) 

while one associated with the second characteristic field is 

H2(C/o): u:^U2{h,Uo)^uo+ ^{h~ho)^lQ^+J^, h>0. (2.11) 

Along the Hugoniot curves 7ii,7i2, the corresponding shock speeds are given by 

hui2 — hoUo 



Ai,2(C/o,C/) = 



h — ho 



As is customary, the shock speed Xi{Uo,U) is required to satisfy Lax shock 
inequalities [17]: 

A,(?7) < A,(C/o,C/) < A,;(C/o), 1-1,2. (2.13) 

Thus, the 1-shock curve Si{Uo) initiating from the left-hand state Uq and consisting 
of all right-hand states U that can be connected to Uo by a Lax shock associated 
with the first characteristic field is 



Si{Uo) ■■ u = ui{h, Uo) = uo- - ho)\j + h > ho. 



(2.14) 



Similarly, the 2-shock curve S2(Uo) issuing from a left-hand state Uo consisting 
of all right-hand states U that can be connected to Uo by a Lax shock associated 
with the second characteristic field is 

52(f/o): u = uiih, Uo) = uo+ ^-{h - ho)^ (1 + -^)' ^ < ^o, (2.15) 
We summarize these results in the following proposition. 
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Lemma 2.2 (Shock wave curves). Given a left-hand state Uq, the l-shock curve 
Si(Uo) consisting of all right-hand states U that can be connected to Uq by a Lax 
shock is 

Si{Uo) ■■ u = ui{h, Uq) = uo - ~ ho)^ + h > Hq. 

The 2-shock curve S2{Uo) consisting of all right-hand states U that can be connected 
to Uq by a Lax shock is 



S2iUo) ■■ u = ui{h, Uo) = uo + - ^o)^ + h < Hq 



In view of the Lax shock inequahties (2.131, we also conclude that the backward 
l-shock curve Si{Uo) issuing from a right-hand state Uq and consisting of all left- 
hand states U that can be connected to Uq by a Lax shock associated with the first 
characteristic field is 



5f(C/o): u^uiih,Uo)^uo- ^{h-ho)^lQ^7J^, h < ho, 



(2.16) 



Similarly, the backward 2-shock curve 5^([/o) issuing from a right-hand state Uo 
and consisting of all left-hand states U that can be connected to Uq by a Lax shock 
associated with the second characteristic field is 

SiiUo): u^u,ih,Uo)^uo + ^ih-~hQ)^lQ^+J^, h>ho, (2.17) 

Next, we discuss the properties of rarefaction waves, i.e., smooth self-similar 



solutions to the system ( 1 . 1 1 associated with one of the two genuinely nonlinear 
characteristic fields. These waves satisfy the ordinary differential equation: 

at VAi • r^{U) 
For waves in the first family, we have 

dh{C) 2h{C) 2 



0. 



It follows that 



du jg 
dh Y /i ' 

therefore, the integral curve passing through a given point Uq — {Hq, uq, oq) is given 
by 

u^Uq- 2y/g{Vh- \/ho). 
Moreover, the characteristic speed should increase through a rarefaction fan, i.e., 

Ai(i7) > Ai(;7o), (2.19) 

which implies 

h> Hq, 
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Thus, we define a rarefaction curve TZi{Uq) issuing from a given left-hand state 
Uq and consisting of all the right-hand states U that can be connected to Uq by a 
rarefaction wave associated with the first characteristic field as 



7^l([/o) : u = vi{h,Uo) ■■= uo -2y/g{Vh~ y^), h < Hq. 
A 1-rarefaction wave is determined by 

2 /X _ Xo\ 

Xt ta) 



Uq 



(2.20) 



(2.21) 



while h is determined by the equation (2.20 1 and the component a remains constant. 

Similarly, we define the rarefaction curve 'R-2{Uq) issuing from a given left-hand 
state J7o and consisting of all the right-hand states U that can be connected to C/q 
by a rarefaction wave associated with the second characteristic field as 

^ (2.22) 



'R2{Uq) : u^V2ih,Uo) ■■=uo + 2^{vh- ^/ho), h>ho 



The w-component of the 2-rarefaction wave is determined by (2.211 and the h- 
component is given by (2.221. 

We can summarize the above results in: 

Lemma 2.3 (Rarefaction wave curves). Given a left-hand state Uq, the l-rarefaction 
curve TZi(Uq) consisting of all right-hand states U that can be connected to Uq by 
a rarefaction wave associated with the first characteristic field is 

TZi{U()) : u = vi{h,Uo) ■.= uo-2^{Vh- \/fH)), h < ho. 

The 2-rarefaction curve Ti-2{Uo) consisting of all right-hand states U that can be 
connected to Uq by a rarefaction wave associated with the second characteristic field 
is 

T^2{U()) ■ u = V2{h,Uo) ■■= Uq + 2^{Vh - \/ho), h>ho. 

We will also need backward curves which we define here for completeness. Given 
a right-hand state Uq, the 1-rarefaction curve TZf{Uo) consisting of all left-hand 
states U that can be connected to Uq by a rarefaction wave associated with the 
first characteristic field is 

^ (2.23) 



7^f (f/o) : u = vi{h,Uo) ■.= uo-2^{Vh^ ^ho), h>ho. 

The 2-rarefaction curve 7?,f ([/q) consisting of all left-hand states U that can be 
connected to Uo by a rarefaction wave associated with the second characteristic 
field is 



7^f (t/o) : u = V2{h,Ua) ■.= uo + 2yqj{Vh^ h < ho. 

In turn, we are in a position to define the wave curves, as follows 
Wl([/o) = 5l(C/o)U7^l(C/o), 
Wf({/o)=5f(^7o)U7^f([/o), 
W2(C/o) = 52(C/o)U7^2(^7o), 

>vf({/o)-52^(^7o)u7^f(^7o). 

Some properties of the wave curves are now checked. 



(2.24) 



(2.25) 
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Lemma 2.4 (Monotonicity properties). The wave curve Wi{Uq) can be parame- 
terized in the form h t-^ u — u{h),h > 0, where the function u is strictly convex 
and strictly decreasing in h. The wave curve W2{U^) can be parameterized in the 
form h ^ u = u{h),h > 0, where the function u is strictly concave and strictly 
decreasing in h. 

Proof. We only give the proof for the 1-wave curve yVi(L'o), the proof for 'W2{Uq) 
being similar. For the shock part iSi(?7o), we have 

1 1 /to 
rfu _ _ [g 2h^ho^ 2fe2 
dh~ V 2 /I 1 ^ ■ 

For the rarefaction part TZi{Uo), we have 

dh \/^ 

This establishes the desired monotonicity property of 'Wi{Uo). 

The convexity of yVi(J7) follows from the fact that du/dh is increasing. Indeed, 
along the shock part 5i(C/o) it holds 



/in \ / 1 I 1 I 1 / 1 I 1 I 'l 



dh? " V 2 



- + — 
h ' ho 



> 



and, along the rarefaction part TZi(Uo), 



— = > n 

dh'^ 2/l3/2 ^ 

which completes the proof. □ 
Next, we consider the 3-curve from a state Uq, which consists of all states U 



that can be connected to Uq by a stationary wave. As seen in (1.9 1, U and Uq are 
related by the Rankine-Hugoniot relations 

[hu] = 

[-+g{h + a)] = 0. 
This leads to a natural definition of a curve parameterized in h: 



(2.26) 



f u^u(h) ^ 

a = a[h) = flQ H + n — tiQ. 

3. Admissibility conditions for stationary waves 

3.1. Two possible stationary jumps. In view of the discussion in the previous 
section, the states across a stationary wave are constraint by the Rankine-Hugoniot 



relations (2.261. From a given left-hand state we have to determine the right 
hand state, which has three components, determined by the two equations ( |2.26 1 



Moreover, since the component a changes only through stationary waves (which 
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propagate with zero speed) for given bottom levels a± we should solve for u and h 
in terms of a. Thus, we rewrite (2.261 in the form 

houo 



ao — a + 



2g 



h — ht 



■0 



0. 



Substituting for u and re-arranging the terms, we obtain 

hoUo 



This leads us to search for roots of the function 

Un / /in 



h — ht 



0- 



Let us set 



.2/1^1/3 



,(L/o) := flo + 



(3.1) 



(3.2) 



1 + /imir 



Some useful properties of the function ip in (3.1 1 are now derived. 

Lemma 3.1. Suppose that Uq = (/io,uo,ao) fl^*^ Q given with uq ^ 0. The 
function ip : (0, +oo) R is smooth and convex and, for some ft-min; H is decreasing 
in the interval (0, /imin) o,nd is increasing in the interval (ft.niin,oo); with 



lim Lp{h) = lim (p{h) ~ +oo. 

h — *0 h — >+oo 



(3.3) 



Furthermore, if a > a,nin then the function ip has two roots h^,{Ua),h* (Uq) with 
h^^Uo) < ft.min(CA)) ^ h*(Uo). Thcsc inequalities are strict whenever a > a,nin(C^o)- 

Proof. The smoothness of the function ip and the limiting conditions are obvious. 
Moreover, we have 



dip{h) 



dh gh^ 
(for Uq =/= 0) which is positive if and only if 



h > 



lhl\V^ 



1 



.(C/o). 



This establishes the monotonicity property of ip. Furthermore, we have 



d^p{h) _ "iulhl 
dh^ ^ gh^ 



> 0, 



which shows the convexity of ip. If a > aniin(C/o)i then iy9(/imin(fA))) < 0- The other 
conclusions follow immediately. □ 



It is straightforward to check: 
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Lemma 3.2. The function /imin satisfies the following inequalities: 



{Uo) > ho, UoeAiU As, 
{Uo) < ho, Uo e A2, 

/imin(t^o) = ho, Uo e C, 



(3.4) 




(3.5) 



ha>h*{Uo) UoeA2. 



(3.6) 



The function amin(t^o) satisfy the following inequalitites: 

Q-niin (C/q) < floj (''-0. M'o) e ^1 u A2 u ^3, 
amin(t/o) = ao, {ho,uo)eC±. 



(3.7) 



The states that can be connected by stationary waves are characterized as follows. 

Proposition 3.3. Fix a left-hand state Uo = {ho,uo,ao) and a right-hand bottom 
level a. 

(i) If Uo and a > aniin(C^o); then there are two distinct right-hand states 



where Ui{Uo) ■= hoUo/hi{Uo),i = 1,2, that can he connected to Uo by a 
stationary wave satisfying the Rankine-Hugoniot relations. 

(ii) If Uo 7^ and a = a^in{Uo), the two states in (i) coincide and we obtain a 
unique stationary wave. 

(iii) Ifuo ^0 and a < amin(C^o); there is no stationary wave from C/_ to a state 
with level a. 

(iv) // Uo = 0, there is only one stationary jump defined by 



We arrive at an important conclusion on stationary jumps. 

Proposition 3.4. For uo ^ 0, the state {hi{Uo),ui{Uo)) belongs to Ai if uq < 0, 

and belongs to A3 if uo > 0, while the state {h2{Uo) , U2{Uo)) always belongs to A2. 
Moreover, we have 



It is interesting to observe that the shock speed in the genuinely nonlinear char- 
acteristic fields will change sign along the shock curves. Therefore, it exchanges its 
order with the linearly degenerate field, as stated in the following theorem. 

Proposition 3.5. (a) IfUo & Ai, then there exists C/q e Si{Uo)r\A2 corresponding 
to h = h> ho such that 



Ui^2 := {hi^2{Uo),ui,2{Uo),a) 



u = Uo = 0, h = ho + a — ao- 




Uq > 0, 
Uo < 0. 



(3.8) 



MUo,Uo)=0, 

~Xi{Uo,U)>0, U €Si{Uo),h€{ho,ho), 
~Xi{Uo,U)<0, U eSi{Uo),he(ho,+oo). 



(3.9) 
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// Uq e A2U A3, then 

~Xi{Uo,U)<0, UeSiiUo). (3.10) 

(b) IfUo e A3, then there exists Uq G (Uq) D A2 corresponding to h — h> Hq 
such that 

~X2iUo,U)>0, U eSiiUo),heiho,ho), (3.11) 
A2(i7o,[/) <0, [/e5f(C/o),/iG (/io,+oo). 
IfUo e AiU A2, then 

MUo,U)>0, UGSiiUo). (3.12) 

3.2. Two-parameter wave sets. From Proposition |3 .4| and the arguments in the 
previous section, we can now construct wave composites. It turns out that two- 
parameter wave sets can be constructed. For definiteness, we now illustrate this 
feature on a particular case. Suppose that Uq = {ho,uo,ao) S Aj. We can use a 
stationary wave from Uq to a state Um = {hm:Um,cim) & ^2" using h* , followed 
by another stationary wave from Um to ?7 e ^1 using the corresponding value h^,, 
then we continue with 1-waves, and as in Ai the characteristic speed is positive. 
As ttm can vary, the set of such states U forms a two-parameter set of composite 
waves containing first and third waves. Such wave sets were constructed even for 
strictly hyperbolic systems by Hayes and LeFloch [12] . 

To make the Riemann problem well-posed, it is necessary to impose an additional 
admissibility criterion. 



3.3. The monotonicity criterion. Since the Riemann problem for (1.1 1 may in 

principle admit up to a one-parameter family of solutions, we now require that the 
Riemann solutions of interest satisfy a monotonicity condition in the component a. 

(MC) (Monotonicity Criterion) Along any stationary curve W3{Uo), the bottom 
level a is a monotone function in h. The total variation of the bottom level 
component of any Riemann solution must not exceed (and, therefore, is 
equal to) \aL — aji\, where aL,oiR are left-hand and right-hand cross-section 
levels. 

A similar selection criterion was used by Isaacson and Temple [TH [TS] and by 
LeFloch and Thanh [21], and by Goatin and LeFloch [Sj. Under the transformation 
(if necessary) 

X — > —X, u — > —u, 

a right-hand state U — {h, u, a) transforms into a left-hand state of the form U' — 
(h, —u,a). Therefore, it is not restrictive to assume that 

flL < flfl- (3.13) 

Lemma 3.6. The Monotonicity Criterion implies that stationary shocks do not 
cross the boundary of strict hyperbolicity. In other words, we have: 

(i) If Uq G U A3, then only the stationary shock based on the value H^^Uq) 
is admissible. 

(ii) If Uq G A2, then only the stationary shock using h*(Uo) is admissible. 
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Proof. Recall that the Rankine-Hugoniot relations associated with the linearly de- 



generate field (2.27) implies that the component a can be expressed as a function 
of h: 

2 2 

a = a[h) = ao H \- h- fiQ, 

2.9 

where 

u = uin) = — , — . 

h 

Thus, differentiating a with respect to h, we find 

uu'{h) houo 
a (hj = h 1 = -u—— + 1 

9 gh^ 

gh 

which is positive (resp. negative) if and only if 

- gh <0 (resp. - gh> Q) 

or {h,u,a) € A2 (resp. £ Ai or G A^). Thus, in order that a' keeps the same sign, 
the point (/i, u, a) must remain on the same side as {ho, uq, ao) with respect to C±. 
The conclusions in (i) and (ii) follow. □ 



It follows from Lemma 3.6 that for a given Uq = (/iqi Wq, ao) G Ai, i = 1,2, 3, and 
a level a, we can define a unique point U = {h, u, a) so that the two points Uq, U 
can be connected by a stationary wave satisfying the (MC) criterion. We have a 
mapping 

SW{., a) : [0,00) xRxR+^ [0, 00) x R x R+ 

Uq = {hQ,uo, ao) i-> SW{Uq, a) = U ^ [h, u, a), ^^-^^^ 

such that Uq and U can be connected by a stationary wave satisfying the (MC) 
condition. Observe that this mapping is single-valued except on the hypersurface 
C, where it has two- values. 

Let us use the following notation: Wi{Uo, U) will stand for the ith wave from a 
left-hand state Uq to the right-hand state U , i = \, 2, 3. To represent the fact that 
the wave Wi{Ui, U2) is followed by the wave Wj{U2, U3), we use the notation: 

W,{Ui,U2)®Wj{U2,U3). (3.15) 

4. The Riemann problem 

In this section we construct the solutions of the Riemann problem, by combining 
Lax shocks, rarefaction waves, and stationary waves satisfying the admissibility 
condition (MC). 

Recall that for general strictly hyperbolic systems of conservation laws, the solu- 
tion to the Riemann problem exist when the initial jump is sufficiently small only. 
That is to say that the right-hand states Ur should lie in a small neighborhood 



of the left-hand state Ul- However, for the system (1.1 1, we can cover large data 
and essentially cover a full domain of existence for any given left-hand state. More 
precisely, we determine the precise range of right-hand states in which the Riemann 
solution exists. 
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4.1. Solutions containing only one wave of each characteristic family. We 

begin by constructing solutions containing only one wave corresponding to each 
characteristic field which is identified as each family of waves. This structure of 
solutions is standard in the theory of strictly hyperbolic system of conservation 
laws. In the next subsection we will consider solutions that contain up to two 
waves in the same family. The following theorem deals with the case where the 
left-hand state Ul is in Ai. 

Theorem 4.1. Let Ul G Ai and set Ui := SW{UL,aii), {U2} = >Vi(C/i) n 
(Ur). Then, the Riemann problem { 1.1 )-{ 1.2) admits an admissible solution 



with the following structure 

WsiUL, Ul) ® WiiUi, U2) e W2{U2, Ur), (4.1) 
provided /12 < ^i- f-Fi^Mrej^. 

Proof. Observe that the set of composite waves S'Vl^(yVi(C/L), a/j) consists of three 
monotone decreasing curves, and each lies entirely in each region Ai,i — 1,2,3. 
The monotone increasing backward curve W;^(?7ij) therefore may cut the three 
composite curves at a unique point, two point, or else does not meet the wave 
composite set. The Riemann problem therefore may admit a unique solution, two 
solutions, or has no solution. 

The state Ul belongs to Ai and in this region, the A3 is the smallest of the 
three characteristic speeds. A stationary wave from Ul = {hL,UL,aL) to Ui = 
{hi,ui,aii) exists, since ql < a^. Moreover, by Lemma [3. 6[ we have Ui € Ai. 

If /12 < hi, then the stationary wave is followed by a 1-rarefaction wave with 
positive speed, and then can be continued by a 2-wave W^2(C^2, C^i?)- If ^2 > hi. 



then the 1-wave in (3.15) is a shock wave. Since ft.2 ^ ^1 and Ui £ Ai, the shock 



speed X2{Ui,U2) > 0, and thus it can follow a stationary wave (with zero speed). 
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Moreover, it is derived from (2.121 that 



Ai(C/i,C/2) = u 



hi 



/12U2 — hiui 
h2 ~ hi 



thus 



Al(f7i,[/2) <U2 + 




This means the 1-shock Si{Ui, U2) can always follow the 2-shock S2{U2, Un)- 
ilar for rarefaction waves. Therefore, the solution structure ( 3.15[ ) holds. 



(4.2) 

Sim- 

□ 



The following theorem deals with the case where the left-hand state Ul is in 

Theorem 4.2. Let Ul G Ai\JA2. Then there exists a region oj values Ur such that 
SW{'Wi{Ul)tO,b) ^y\^2 {Ur) 7^ ^- ^his case this intersection may contain either 
only one or both points Ui € A2 and U2 G ^3. The Riemann problem (1.1)-1.2 
therefore has a solution with the structure 

Wi{Ul, U3) ® WsiUs, Ul) ® W2{Ui, Ur), (4.3) 

where U3 £ Wi(C/l) is the point such that Ui — SW{U3,aR), and also 

Wi{UL,Ui) e WsiUi, U2) e W2iU2, Ur), (4.4) 

where U4 € Wi(Ul) is the point such that U2 ~ SW{U/i,aR), if h2 > hR whenever 
UrG A^. (Figure 3 ) 



Proof. The solution may begin with a 1-wave, either 1-shock with a negative shock 
speed to a state U^, or a 1-rarefaction wave with \i{Uz) < 0, followed by a station- 
ary wave W3{U3,Ui) from C/3 to Ui, then followed by a 2-wave W2{Ui,Ur) from 
Ul to Ur. It is similar in the case of J72. However, in order that the stationary 
wave W3{U4,U2), for some U4 S yVi{UL) and U4 G A^ obviously, to be followed 
by a 2-wave W2{U2, Ur), it is required that the wave is a shock with non-negative 
shock speed X2{U2, Ur). This is equivalent to h2 > Tir. □ 

Theorem 4.3. Let Ul G A3 and Ur G A1UA2, and set Ui = SW{'Wi{UR),aL)n 
Wi{Ul) and U2 = SW{Ui,aR,) G WHUr). 



(i) //C/i G ^2 UC+ U {m = 0}, the Riemann problem (l.l)-(l. 2) has a solution 
with the following structure 

Wi{Ul,Ui) e W3{Ui, U2) e W2iU2, Ur). (4.5) 



(ii) Lf Ul G A.2 UC_, provided Hr > h2, the Riemann solution (4-5) also exists 



(iii) Lf Ul G Ai U A3, the construction (4.5) does not make sense. 
(Figure^ . 
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Figure 3. Solution for Ul £ AiU A2 




Figure 4. Solution for Ul E A3 



Proof. If Ul E A2U C, the non-positive speed wave Wi{Ul, Ui) can be followed by 
a stationary wave W3{Ui, U2). 

When U2 G A2 — , if C/i € A2 LI C+ U {u — 0}, then this stationary wave can 
always be followed by a 2- wave W2{U2, Uji), since the wave speed of the 2- wave is 
positive. This establishes (i). 

If Ul G A2 UC_, the wave speed of the 2-wave W2{U2,Ub) is non-negative if 
and only if hR>h2- This proves (it). 
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Figure 5. Ul may be anywhere 



If Ui G Ai, the 1-wave has positive speed. So it can not be followed by a 
stationary wave. If Ui S A3, then U2 G A3 by the (MC) criterion. So the 2-wave 
W2{U2, Ufi) has negative speed. So it can not be proceeded by a stationary wave. 
This proves (iii). □ 

The above theorem enables Ur to vary in each region Ai,A2, and A3. The next 
theorem enables Ul to vary in all the three regions. 

Theorem 4.4. Let Ur S A3. Set Ui = SW{UR,aL), U2 = Wi{Ui) n Wi{Ul). A 
Riemann solution exists and has the following structure 

Wi{Ul,U2) ® W2{U2, Ul) ® W3{Ui, Ur), (4.6) 
provided h2 <hi. (Figure^. 

Proof. The stationary wave W3{Ui, Ur) turns out to have the greatest wave speed. 
In order for this wave to be proceeded by the 2-wave W2{Ui,U2), the wave speed 
of this 2-wave has to be non-positive. This is equivalent to the condition h2 < hi. 



according to Proposition (3.5). Similar to (4.4), we have 



Al(f/L,i72) < A2(C/2,C/l). 

so that the 1-wave Wi{Ul, U2) can follow the 2-wave W2{U2, Ui). □ 

4.2. Solutions containing more than one wave of each characteristic fam- 
ily. It is remarkable feature of the shallow water system that we can also construct 
solutions with four elementary waves, using three available characteristic fields. 
This illustrates one of the difficulties in coping with the Riemann problem when 
the system under consideration is not strictly hyperbolic. 
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Figure 6. Solution with repeated two 1-waves 



Theorem 4.5. Let Ul e A2UA3 and set U+ = Wi{UL)nC+, Wi] ^ SW{U+,aR)n 
Ai,{U2} = yVi{Ui) n WHUr). The Riemann problem ^7^-{1.2) has a solution 
with the following structure 

RiiUL, U+) ® 1¥3(C/+, f/i) © Wi{Ui, U2) ® W2{U2, Ur), (4.7) 

provided /12 < ^i- ffi^Mrej^. 

Theorem 4.6. For any Ul, set {Ui} = SW{C-,aR) n WHUr) n A2, U2 = 
{h2,U2,aL) e C_ such that Ui ^ SW{U2), and {U3} = Wi{U2) n Wi{Ul). Then 
the Riemann problem (1.1)-(1.2) has a solution with the following structure 

Wi{Ul, U3) e i?2(C/3, C/2) ® W3iU2,Ui) © W2iUi, Ur), (4.8) 

provided Hr > hi and < /i2. fFi(7Mre[^. 

Thus, we see that the Riemann problem ( |l.l[ )- pr2| has a solution consisting of 
a 1-, a 3-, and two 2-waves. 

It is interesting to note that there are solutions satisfying the (MC) criterion 
which contain three waves with the same speed (zero). This is the case when a 
stationary wave jumps from the level a = to an intermediate level am between 
aL and aR, followed by an "intermediate" fc-shock with zero speed at the level a™, 
k — 1,3, and then followed by another stationary wave jumping from the level a„ 
to gr. Thus, there are only two possibilities: 

(i) Ul belongs to Ai and a 1-shock with zero speed is used. 

(ii) Ur belongs to A3 and a 2-shock with zero speed is used. 

We just describe the first case (i), as the second case is similar. Recall from 
Proposition 3.5 that for any U €z Ai, there exists a unique point denoted by C/ G 
Wi{U) n A2 such that 

Ai([/,C/) =0. 
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Figure 7. Solution with repeated two 2-waves 

Theorem 4.7. Let Ul G Ai and set 

SWiUL, [aL,aR]) := U,e[„,,„,]^W^([/L, a), 
SWiULdaL,aR]) :={C/| f/ G 5W^(C/l, [a^, a^])}. 

Whenever 

SWiUL, [aL,aR]) n WHUr) + 
there exist a„j G [flLjafl]; V\ = ^^^(C/i, a„j), anrf 

that defines a solution with the structure 

W3{Ul, Ul) © Si{Ui,Ui) © W^iUi, U2) © W2{U2, Ur). (4.9) 
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